A lattice study of the strangeness content of the nucleon 
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Abstract 

We determine the quark contributions to the nucleon spin As, Au and Ad as well as their 
contributions to the nucleon mass, the cr-terms. This is done by computing both, the quark 
line connected and disconnected contributions to the respective matrix elements, using the non- 
perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate np = 2 
mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a ~ 0.073 fm. 
The renormalization of the matrix elements involves mixing between contributions from different 
quark flavours. The pion- nucleon cr-term is extrapolated to physical quark masses exploiting the 
sea quark mass dependence of the nucleon mass. We obtain the renormalized value ct^n = (38 ± 12) 
MeV at the physical point and the strangeness fraction f Ts = a s /m N = 0.012(14)1^° at our larger 
than physical sea quark mass. For the strangeness contribution to the nucleon spin we obtain 
As^(\/7T4GeV) = -0.020(10)(2). 

1 Introduction 

Most of the nucleon's mass is generated by the spontaneous breaking of chiral symmetry and only a 
small part can be attributed directly to the masses of its valence and sea quarks. The quantities 

h q = m g (N\qq\N)/m N = a q /m N (1) 
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parameterize the fractions of the nucleon mass m^ that are carried by quarks of flavour q. These scalar 
matrix elements also determine the coupling strength of the Standard Model Higgs boson (or of any 
similar scalar particle) at zero recoil to the nucleon. This then might couple to heavy particles, some of 
which are dark matter candidates pQ. The combination ran fr q , q £ {u, d, s}, appears quadratically 
in the cross section that is proportional to |/n| 2 , where 



m N 



5>.— 

Z — ' Vfln 

q£{u,d,s} 



9n, 



(2) 



q£{c,b,t,. 



with the couplings a q oc m q /mw n h denotes the number of heavy quark flavours. Note that due to 
the trace anomaly of the energy momentum tensor one obtains 



fr q , 
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so that the coupling /n only depends mildly on the number and properties of (discovered and not yet 
discovered) heavy quark flavours [2]. 

The light quark contribution, the pion-nucleon cr-term, is defined as 
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From dispersive analyses of pion-nucleon scattering data, the values [3] g^n = 45(8) MeV and [1] 
cr^ = 64(7) MeV were obtained while a recent covariant baryon chiral perturbation theory (B%PT) 
analysis of the available scattering and pionic atom data [5] resulted in the estimate ovn = 59(7) MeV. 

Not only do the quarks contribute a tiny fraction to the nucleon's mass, even the nucleon spin 
is mostly carried by the gluons. This spin can be factorized into a quark spin contribution AS, a 
quark angular momentum contribution L q and gluonic contributions AG and Lq (for spin and angular 
momentum): 

^ = ^AS + L g + AG + L G . (5) 

In the naive non-relativistic SU(6) quark model, AS = 1, with vanishing angular momentum and 
gluon contributions. In this case there will be no strangeness contribution As to AS = Au + Ad + 
As + ■ ■ • where, in our notation, Ag contains both, the spin of the quarks q and of the antiquarks 
q. Experimentally, As is obtained by integrating the strangeness contribution to the spin structure 
function g\ over the momentum fraction x. The integral over the range in which data exists usually 
agrees with zero, see e.g. new COMPASS data [6] for x > 0.004, while global analyses tend to obtain 
values As w -0.12 

Here, we directly compute the quark line disconnected (and connected) contributions to the scalar 
and axial matrix elements that appear in the above quantities [9j [10]. Other recent direct lattice 
determinations of these quantities include Refs. [HI [121 EES EEH EEIS [16]. We will first describe the 
methods used, before we present results on the quark spin contributions and the scalar matrix elements. 



2 Methods and simulation parameters 

We simulate n-p = 2 non-perturbatively improved Sheikholeslami-Wohlert Fermions, using the Wilson 
gauge action, at (3 = 5.29 and k = n u d = 0.13632. Setting the scale from the chirally extrapolated 
nucleon mass [17], we obtain the lattice spacing a -1 = 2.71 (2) (7) GeV, where the errors are statistical 
and from the extrapolation, respectively. 
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Figure 1: Dependence of R dls on tf for the scalar strangeness matrix element on the 40 3 64 volume for 
smeared-smeared (SS) and smeared-point (SP) two-point functions, together with the fit result. 



We realize two additional valence k values, n m = 0.13609 and k s = 0.13550. The corresponding pseu- 
doscalar masses are m PSM = 285(3)(7) MeV, m PS , m = 449(3)(11) MeV and m PS , s = 720(5) (18) MeV. 
The strange quark mass was fixed so that the mps, s value is close to the mass of a hypothetical strange- 
antistrange pseudoscalar meson: (m 2 K± +m 2 K0 — m^ ± ) 1//2 w 686.9 MeV. We investigate volumes of 32 3 64 
and 40 3 64 lattice points, i.e., Lmps — 3.36 and 4.20, respectively, where the largest spatial lattice extent 
is L « 2.91 fm. 

The matrix elements of interest are extracted from the large time behaviour of ratios of three-point 
functions over two-point functions where we create a proton at a time t — and destroy it at zero 
momentum at a time tf. At an intermediate time < t < tf the current of interest is inserted. The 
three-point function contains both, a quark line connected and a quark line disconnected contribution. 
For the example of the axial current the ratios to be calculated read 

.con. f) = (rgQf,*)) dl = (rg 1 <"(^)ExTr(7 J 7 5 M-Hx,t;x,t))) 

\ L unpol L '2ptl t f )l \ l unpol U 2pil t f )) 

where M is the lattice Dirac operator, r unpo i = (1 + 74)/2 a parity projector and r po i = ^7j75r unpol 
projects out the difference between the two polarizations (in direction j). We average over j = 1,2,3 
to increase statistics. For the scalar case we have to replace 7^75 >-)■ 1, r po i i-> r unpo i and add the 
vacuum condensate Tr M _1 (x, t; x, t)) to R dls above. For the up and down quark matrix elements 
we compute the sum of connected and disconnected terms while only R dls contributes to As and a s . 

The disconnected contribution is computed with the methods described in [181 El EE] where we fix 
t = 4a ~ 0.29 fm and vary tf. Employing optimized sink and source smearing we find the asymptotic 
limit to be effectively reached for tf > 5a and fit the ratios to a constant for tf > 6a w 0.44 fm, see 
Fig. [T]for an example. 



3 The quark contributions to the proton spin 

Non-singlet axial currents renormalize with a renormalization factor Z^{a) that only depends on the 
lattice spacing. This was determined non-perturbatively for the action and lattice spacing in use [i~9] : 
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Figure 2: Volume and (light) valence quark mass dependence of the unrenormalized As lat . 

Table 1: The connected and disconnected contributions to Ag lat of the proton on the 40 3 64 
volume (L 2.91 fm) as well as the renormalized spin content at a scale « ~ V7.4 GeV. 
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0.76485(64) (73). However, due to the axial anomaly, the renormalization constant of singlet 
currents, Z A (fi, a), acquires an anomalous dimension [20J and will depend on the scheme and scale used 
in the continuum. This has been determined perturbatively [2T] and the result for the conversion to 



the MS scheme reads 

z(fi,a) := Z s A (fi,a) - Z n A s [a) = C F n F [15.8380(8) - 61n(aV 



4tJ 



+ 0(a 3 s ) . (7) 
1 and allow for a 50 % 



We extract an improved coupling from the measured plaquette, set ft = a 
systematic error on 2(a/7.4 GeV) = 0.0055(1)(27). As can be seen from the small anomalous dimension, 
the scale dependence of z(fi) is quite mild. Perturbative 0(a) improvement is implemented [22] to 
0(a), where again we allow for a 50 % systematic error. 

In the n-p = 1 + 1 + 1 theory the matrix elements renormalize as follows. 

AY^ifi) = (Am + Ad + As^ifi) = Z s A {fi, a)(Au + Ad + As) lat (a) , (8) 
a 8 = AT 8 = {Au + Ad- 2As) ME = Zl s (a){Au + Ad - 2As) lat (a) , (9) 
g A = AT 3 = (Au - Adf^ = Z n A s (a)(Au - Ad) lat {a) . (10) 



We employ n-p = 2 mass-degenerate sea quarks so that our singlet current is Am + Ad instead. This 
modifies the renormalization pattern [TO] 
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where As MS receives light quark contributions but the Au MS and Ad M remain unaffected by the 
(quenched) strange quark. We remark that unitarity is violated, due to the partial quenching. 

In Fig. [2] we display the volume and (light) valence quark mass dependence of our unrenormalized 
As lat . There are no significant finite size effects and we take the independence on the valence quark 
mass as an indication that our result may also approximately apply to physical light quark masses. We 
display our 0(a) improved results with statistical and systematic errors in Table [TJ The Au MS and 
Ad MS values are reduced by about 0.035, due to the sea quark contributions while As MS increases by 
0.002 (< 10 %), due to the mixing with light quark flavours. 

Note that we find Qa ~ 1.11(2). This underestimation of the value g& = 1.267(4) from neutron 
/3-decays can probably be explained by our twice as heavy as physical pseudoscalar meson. Our main 
finding is a small negative As ms (a/7.4 GeV) = —0.020(11) that is unlikely to decrease significantly 
if the sea quark mass is reduced: the mixing effects on the renormalization are small, in spite of the 
comparatively large Am and Ad values, and so is the dependence on the light valence quark mass, see 
Fig. El 



4 The light and strange cr-terms 

In massless schemes that preserve the chiral symmetry the combinations m q qq are invariant under 
scale transformations (up to lattice artefacts). This operator identity holds for expectation values, 
independent of the external state. The Wilson action however explicitly breaks the chiral symmetry 
so that singlet and non-singlet flavour combinations will renormalize differently. It turns out that at 
our lattice spacing the ratio between the two renormalization factors deviates from unity by as much 
as 40 %. Consequently, the renormalized strangeness matrix element receives large subtractions from 
light quark contributions. For details, see [H] and references therein. Moreover, without considering 
mixing with the gluonic operator aGG (that we have not determined), we are not able to implement 
O(a) improvement. Note that due to the trace anomaly of the energy-momentum tensor such gluonic 
operators also become relevant for heavy quark masses in the continuum theory [2]. We plan to take 
this effect into account in the future. 

Again we find no significant finite size effects and obtain upsn(^ps ~ 285 MeV) = 106(H)(3) MeV. 
Recently, the ETM Collaboration [TTJ reported the value a PSN (m P s « 380 MeV) = 150(1)(10) MeV, 
using a different lattice action. If we assume the leading order chiral behaviour dm^/dmp s = const then 
Eq. (j3J) suggests to rescale the results according to the ratio of the squared pseudoscalar masses ~ 1.78. 
However, 1.42(19) < 1.78: higher order chiral corrections are relevant. Using additional nucleon mass 
data we extrapolate our value to the physical point [T7] and obtain [5] 

a^ ys = (38±12)MeV, (12) 

where the dominant error is from the chiral extrapolation, see Fig. [3J 

We are also able to compute the strangeness and gluon contributions to the nucleon mass, see 
Eqs. (pQ) and ©: 

f Ts = 0.012(14)^° , f Ta = 0.95111?. (13) 

The light and strange quarks contribute a fraction between 3 % and 8 % to the nucleon mass. The large 
uncertainty on fx s or, equivalently, a s = 121^ MeV is due to large cancellations in the renormaliza- 
tion. To reduce these, with Wilson Fermions, one will need to simulate at finer lattice spacings where 
differences between singlet and non-singlet renormalization constants become smaller. We are able to 
state a 95 % confidence level upper limit on the phenomenologically relevant y-ratio: y < 0.14. 
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Figure 3: Extrapolation of er PSN /mp S to the physical point [T7] using covariant B%PT for the 40 3 64 
volume (solid symbol). The broad error band is obtained from nucleon mass data alone. The horizontal 
line is the leading order expectation and the open symbol our result for the 32 3 64 volume. 

5 Conclusion 

We calculated disconnected contributions to the proton structure. We find a pion-nucleon cx-term 
ct^tv = 38(12) MeV at the physical point and an upper limit on the ratio of scalar strangeness over light 
quark content of y < 0.14. The renormalized light sea quark contributions amount to less than 10 % of 
the total matrix elements, both for the spin content and for the cr-term. However, at our lattice spacing, 
prior to the renormalization, these account for 30 % of the bare scalar lattice matrix elements [9] and 
need to be taken into account for Wilson Fermions. 



Our small negative value As MS (\/7A GeV) = —0.020(11) indicates stronger than naively expected 

violations of SU(3).f symmetry in weak decays. This impacts on determinations of polarized parton 
distribution functions [231 13 IB] where the constraint on the integrals by the (assumed) proton tensor 
charge value a 8 = 3F — D should probably be relaxed. A small value of As was also reported in [13], 
albeit without renormalization. In view of these findings, lattice studies of the spin content of hyperons 
and of their weak transition matrix elements seem particularly interesting. 
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